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Abstract In this brief summary I describe our recent 
work on superfluid neutron star dynamics. I review re- 
sults on shear viscosity, hyperon bulk viscosity, vortex 
mediated mutual friction and the modelling of multifluid 
systems in general. For each problem I provide a set of 
questions that need to be addressed by future work. 
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1 Introduction 

As the temperature of a system decreases towards ab- 
solute zero, matter either freezes to a solid or becomes 
superfluid. The first case is dominated by a structured 
lattice. In the second case, the system acts as a macro- 
scopic quantum system. Regardless of the outcome, the 
rules that apply in an extremely cold system may be very 
different from those that describe the system at higher 
temperatures. This is, of course, a well-known fact and 
the different possible phases of matter have been studied 
in great detail. In particular, increasingly sophisticated 
laboratory experiments continue to provide insights into 
the details of superfluid/superconducting systems, rang- 
ing from the standard Bose-Einstein condensates to su- 
perfluid Helium (Donnelly 1991) and atomic systems ex- 
hibiting fermion Cooper pairing. This area of research is 
highly relevant for those that are fascinated by neutron 
star physics. In fact, mature neutron stars may provide 
the ultimate testing ground for theoretical physics. They 
are expected to contain an elastic nuclear crust perme- 
ated by superfluid neutrons. An outer core where su- 
perfluid neutrons coexist with superconducting protons 
transitions to the deeper core which may contain exot- 
ica like superfluid hyperons and deconfined quarks in a 
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colour superconducting state. As if this was not enough, 
the presence of a "solid" core remains a possibility (even 
though the relevant lattice parameters are unknown). 

Trying to understand the details of the various neu- 
tron star phases and their possible impact on astronom- 
ical observations is an exciting challenge. First of all, it 
requires a working knowledge of much of modern the- 
oretical physics. Secondly, one will need some under- 
standing of the observations. This includes decades of 
data for radio pulsars, in particular associated with the 
glitch events that provide the strongest evidence of the 
presence of (at least partially) decoupled superfluid com- 
ponents. X-ray observations of accreting neutron stars 
in binary systems, and isolated neutron stars closer to 
the Earth, are common and the data for the existence 
of magnetars is now very convincing. In fact, the re- 
cent observations of likely crustal oscillations following 
the giant flares in SGR1806-20 and SGR1900+14 [see 
the contributions by Israel, Watts and Samuelsson] may 
be the harbingers of neutron star asteroseismology. Fu- 
ture gravitational-wave observations should add another 
observational window, although it will likely require an 
advanced generation of detectors sensitive at high fre- 
quencies. 

As the observations continue to improve we will have 
more precise quantitative opportunities to test our the- 
oretical models. To meet these tests the current models 
have to be improved. In fact, many of our "favourite" 
models are far too rudimentary to pass any closer scrutiny. 
This is not surprising given that the relevant theory 
problems are difhcult. Consider for example the possibil- 
ity of superfluid vortex pinning in the neutron star crust. 
At the phenomenological level, the association between 
the observed spin-up of the crust in a glitch and the 
transfer of angular momentum from an initially faster 
spinning (superfluid) component is natural. To invoke 
vortex unpinning as the key agent that motivates the 
event is also natural. Yet it is clear from the effort that 
has gone into glitch modelling that it is difhcult to de- 
scribe the process in a truly quantitative way. 
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The dynamics of a superfluid neutron star is asso- 
ciated with a number of similar problems that need to 
be better understood. The main aim of the research de- 
scribed in this brief review is to arrive at such an un- 
derstanding. Right now we are quite far away from this 
goal. In some cases, e.g. problems concerning basic shear 
viscosity, we may have a clear picture already, but in 
others, eg. how the calculated shear viscosity coefficients 
are used in the case of a superfluid system, we are only 
beginning to be able to pose (what might be) the right 
questions. 

This brief review is a summary of my presentation 
at the London neutron star meeting. It is in no sense a 
complete review of the various problems. Rather it is a 
self-centered view of my thinking at the time of writing. 
The interested reader will find references to the relevant 
literature in the cited papers. 

Before I discuss our recent work, let me provide a list 
of problems that I find particularly interesting: 

How do pulsar glitches really work? As already al- 
luded to above, this remains a vexing issue. In my (some- 
what pessimistic view) view, the present models are not 
truly quantitative. They likely contain the right elements, 
but the details of the underlying physics are not yet 
agreed upon. It is also worth noting that most models 
predict how the glitch proceeds and how the system re- 
laxes back to quasi-cquilibrium. Very few models provide 
a real mechanism to explain why the glitch happened in 
the first place. 

How do we understand neutron star free precession? 
The observational indications that some neutron stars 
are wobbling now seem relatively strong. The obvious 
question is why this behaviour is so rare. After all, free 
precession is the most natural mode of motion for any 
rotating body. The answer will provide insight into the 
damping of fluid motion in the star, and thus may con- 
strain our models for internal dissipation. There could 
also be implications for vortex pinning or vortex-flux 
tube interactions in the superfluid/superconducting core 
[see Link's contribution]. It would also be nice to have 
a model for the excitation of the precession in the first 
place. 

Do the superfluid degrees of freedom affect the oscil- 
lation properties of the star? This is an important ques- 
tion that impacts on future attempts to probe neutron 
star physics via asteroseismology. To provide a "useful" 
answer we need to study the problem using realistic mod- 
els for the superfluid pairing gaps (which, of course, are 
not agreed upon by theorists) within general relativity. 
The latter is key in order to make the results accurate 
enough that it is meaningful to compare to observations. 
In principle, I don't think the mode calculation would be 
very difficult. At least not for non-rotating stars, and as 
long as one is prepared to accept that the pairing gaps 
are likely to remain unknown up to perhaps a factor of 
a few. However, if one wants to account for the presence 
of the crust (which should be penetrated by superfiuid 



neutrons) and the magnetic field, then the modelling be- 
comes miich more challenging. In the case of the crust, 
we have not (until very recently (Carter & Samuelsson 
2006)) had any theoretical formulations that allow for 
the presence of a superfluid component, while in the mag- 
netic field case we do not have a good understanding of 
the internal field structure. 

Neutron stars may radiate detectable gravitational 
waves through a number of scenarios. One possibility is 
that the inertial r-modes are driven unstable as they radi- 
ate gravitationally. The basic mechanism behind this in- 
stability is well understood, and a large number of damp- 
ing mechanisms that counteract the instability have been 
considered. Thus it has become clear that the key decid- 
ing factors relate to superfluidity. The presence of hy- 
perons in the deep stellar core may prevent the instabil- 
ity from happening, but if the hyperons are superfiuid 
then the chemical reactions that lead to bulk viscosity 
are suppressed and the effect on the instability may not 
be considerable. Another important mechanism is the so- 
called mutual friction, which in this context relates to the 
damping due to the scattering of electrons off of magnetic 
fields associated with the superfluid neutron vortices. So 
far all studies of mode-damping due to this effect have 
been based on straight vortices (as in a rigidly rotat- 
ing star). But this may not be the appropriate model! 
If there is a significant oscillation in the star, which has 
some vorticity associated with it, then it seems likely 
that the vortices will get tangled up. The system would 
then be in a "turbulent" state, and from the analogous 
problem in superfluid Helium (Donnelly 1991) we know 
that the mutual friction force is then rather different. 

Finally, I would like to emphasise the multi-fluid as- 
pects of these problems. In most available studies, con- 
cerning for example viscosity, it is implicitly assumed 
that we can apply the standard single fluid equations 
of motion. Yet a multi-fluid system has extra degrees 
of freedom (the second sound in Helium and the analo- 
gous "superfluid modes" in an oscillating neutron star). 
These arc unlikely to be "passive" . In fact, we know that 
the equations that are used to model simple superfluid 
neutron stars admit a so-called two-stream instability 
(Andersson & Comer 2004) . One can speculate that this 
instability becomes relevant when the two components of 
the star rotate at different rates, as in the case of a spin- 
ning down neutron star with a pinned superfluid com- 
ponent. The instability mechanism is very simple and 
familiar from othc;r physical systems, and it would be 
interesting to know whether it can operate in a neutron 
star as well. 



2 Shear viscosity 

We have recently investigated the effects of superfluidity 
on the shear viscosity in a neutron star core (Andersson, 
Comer & Glampedakis 2005). We were motivated to do 
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this by a puzzling result from the literature. The avail- 
able results suggest that the shear viscosity is stronger in 
a superfluid neutron star than it is when the star forms 
a normal-fluid system (Cutler & Lindblom 1987). This 
result contradicts our experience from other superfluid 
systems like Helium, and since the shear viscosity af- 
fects mode damping, free precession and relaxation after 
spin-up events in a crucial way, it is important that we 
understand it. 

To address the problem we combined existing theo- 
retical results for the viscosity coefficients with data for 
the various superfluid energy gaps into a "complete" de- 
scription. This lead to a simple model for the electron 
viscosity which is relevant both when the protons form 
a normal fluid and when they become superconducting. 
This turns out to be the key distinguishing factor. 

Below the neutron superfluid transition temperature 
the dominant contribution to the shear viscosity comes 
from the scattering of relativistic electrons. To model 
the shear viscosity we need to consider two contribu- 
tions. Above the transition temperature at which the 
protons become superconducting, electron-proton scat- 
tering leads to a viscosity coefficient 



(1) 



where Xp is the proton (electron) fraction, Ts — T/IO^K 
and pi5 = p/lO^^g/cm'^. Meanwhile, when the protons 
are superconducting the dominant effect is due to elec- 
trons scattering off of each other. Then we have 



(2) 



The protons play the key role since individual scattering 
processes add like "parallel resistors" . That is, we have 
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where 



» Te, 



(3) 



and it is clear that the most important contribution 
comes from the most frequent scattering process. 

When the protons are superconducting the electron- 
proton scattering will be suppressed, essentially because 
there will be fewer states available for the protons to 
scatter into. In order to allow for the transition to proton 
superconductivity, we can introduce a suppression factor 



TZp such that 



' cp 



' ep 
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Far below the critical transition temperature at which 
the protons become superconducting we should have TZp — 
0, and we see from (3) that the electron-electron scatter- 
ing then dominates the shear viscosity. 

Our results, which are illustrated in Figure 1, explain 
in a clear way why proton superconductivity leads to a 
significant strengthening of the shear viscosity. It should 
be noted that the superfluidity of the neutrons is not 
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Fig. 1 Various shear viscosity estimates for a neutron star 
at temperature 10* K. The solid blue line shows our estimate 
r^ee while the dashed blue curve corresponds to r/ep. These 
estimates are compared to i) the total shear viscosity in a 
normal fluid neutron star, which is dominated by neutron- 
neutron scattering (FI79/CL87, solid red line) taken from 
Cutler & Lindblom (1987), and ii) results obtained by Flowers 
and Itoh for the electron-electron shear viscosity in the crust 
region (crust, FI76, dotted blue curve). The crust is indicated 
by the grey region. [See Andersson, Comer & Glampedakis 
(2005) for the relevant references.] 



the factor which leads to electron-electron scattering be- 
coming the main shear viscosity agent. Rather, it is the 
fact that the onset of superconductivity suppresses the 
electron-proton scattering. As can be seen from Figure 1, 
the electron-electron shear viscosity is not too different 
from the result for neutrons scattering off of each other. 
This means that, in the temperature range where shear 
viscosity dominates, the damping of neutron star oscil- 
lations will be quite similar (modulo multi-fluid effects) 
in the extreme cases when i) the neutrons and protons 
are both normal, and ii) when the neutrons and protons 
are superfluid/superconducting, respectively. The con- 
trast with the case when the neutrons are superfluid and 
the protons normal (and viscosity is dominated by Tyop) 
is clear from Figure 1. This is an interesting observa- 
tion because it shows that proton superconductivity (or 
rather absence thereof) could have a signiflcant effect on 
the dynamics of a neutron star core. 

Future work needs to i) implement the different de- 
grees of freedom of a multifiuid system, and ii) determine 
the required suppression factors (like TZp). 



3 Hyperon viscosity 

The presence of hyperons is expected to affect a neu- 
tron star in several important ways. First of all, the 
carries negative charge which means that the lepton frac- 
tions drop significantly following its appearance. In fact, 
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Fig. 2 Various estimates of the r-mode instability for a sim- 
ple neutron star model with a hyperon core. The figure shows 
the critical rotation rate O/Hk as a function of the core tem- 
perature T. The solid blue curve represents the case of normal 
fluid hyperons. The r-modes are unstable above the curve, 
and even though there appears to be a sizeable instability 
window the hyperon viscosity provides a strong suppression 
(compare to the solid black curve which represents the case 
when there are no hyperons in the star). Basically, the re- 
maining instability window is at such high temperatures that 
the star would cool through the instability regime in a minute 
or so. Far too fast to allow the modes to grow to a large am- 
plitude. The dashed blue curve accounts for rotational effects, 
and shows that the modes are unstable in the fastest spin- 
ning stars also at much lower temperatures. The solid red 
and dashed black curves assume that the hyperon viscosity is 
suppressed below a given critical temperature (here taken to 
be low, at 10^ K), and correspond to the case when there is 
a crust (and a viscous Ekman layer is in operation) and the 
star is entirely fluid, respectively. 



in some models there are virtually no electrons present 
in the core of the star. This may have implications for 
the shear viscosity results that we discussed earlier. Sec- 
ondly, the hyperons can act as an extremely efficient re- 
frigerant. This is mainly because the hyperons may un- 
dergo direct URCA reactions essentially as soon as they 
appear, in clear contrast to the protons which must ex- 
ceed a threshold value of Xp ~ 0.1. Studies of the effect 
of this enhanced cooling have shown that a neutron star 
with a sizeable hyperon core would cool extremely fast. 
In fact, it would cool so fast that it would be much colder 
than observational data suggests. This could be taken as 
evidence against an exotic core, but more likely it sug- 
gests that the hyperons are (at least partly) superfluid. 
This is a natural explanation since superfluidity leads to 
a reduction of the relevant reaction rates. 

Analogous non-leptonic reactions have recently been 
invoked as a mechanism for producing a very strong 
bulk viscosity that may suppress the gravitational-wave 
driven instability of the r-modes (Lindblom & Owen 2002) 
[see also the contribution by Chatterjee & Bandyopad- 
hyay] . This mechanism is particularly important because, 



in contrast to the standard bulk viscosity associated with 
/3-reactions it is relevant at low temperatures (the co- 
efficient scales as rather than T^). For a while it 
was argued that the presence of hyperons would make 
the r-mode instability completely irrelevant. This is now 
understood not to be the case (Nayyar& Owen 2006). 
Basically, the hyperons must be superfluid in order not 
to lead to conflicts with cooling data. Then the nuclear 
reactions that lead to the bulk viscosity must also be 
suppressed, and the effect on the r-mode instability may 
not be so great after all, see Figure 2. In fact, the most re- 
cent estimates suggest that the r-modes may be unstable 
in accreting neutron stars in LMXBs, possibly regulating 
the spin of these stars. Since the associated gravitational- 
wave signal may be detectable by advanced detectors, it 
is important that we understand this problem better. 

Future work needs to i) improve our understanding 
of the hyperon pairing gaps and the associated suppres- 
sion factors for bulk viscosity, and ii) develop a multi- 
fluid description of bulk viscosity that can be used below 
the relevant superfluid transition temperature. 



4 Modelling multifluids 

The equations that govern a general conservative mul- 
tifluid system can be derived from a constrained varia- 
tional principle (Prix 2004; Andersson & Comer 2006a). 
The fundamental variables in this framework are the 
number densities rix, where x is a "constituent" index 
that identifies the different particle species, and the asso- 
ciated transport velocities v]^. To complete the model one 
must provide an energy functional which represents the 
equation of state. In the isotropic (meaning that there 
are no prefered directions) two-fluid problem this func- 
tional takes the form E{nn, rip, w^) where wj^ = vj — vf 
( X ^ y are constituent indices). This immediately leads 
to 
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(5) 



(6) 



(with me the baryon mass) defines the chemical poten- 
tial per unit mass jl^ and the entrainment parameter a. 

The equations of motion for the system can be writ- 
ten 



dtu^ + V^in^vl) = 
and 



/^x - 2™BW> 



= 



(7) 



(8) 



where represents the Lie derivative along (see 
Andersson & Comer (2006b) for an explanation why the 
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Lie derivative is natural in this context). In (8) the mo- 
mentum (per particle) is given by — Ta^{v^ + Ey^w^'), 
where Ex = '^ajm^ny^. From this we can understand the 
nature of the entrainment effect. It is such that the veloc- 
ity and momentum of each constituent in the multifluid 
system are no longer parallel. An alternative, perhaps 
more intuitive, description would represent this mecha- 
nism in terms of an effective mass m* (Andersson, Sidery 
& Comer 2006). In the case of the outer core of a neu- 
tron star, where superfluid neutrons coexist with super- 
conducting protons the entrainment arises because of the 
strong interaction. Each proton/neutron is endowed with 
a cloud of particles of the opposite species, and when it 
flows it drags some of the other fluid along with it thus 
affecting its momentum. 

The entrainment is known to play an important role 
in neutron star dynamics. In particular, it has been shown 
to affect the spectrum of stellar pulsation (Prix, Comer 
& Andersson 2004; Gusakov & Andersson 2006). More- 
over, it is one of the key ingredients in the prescription 
for the vortex mediated mutual friction. 

Future work needs to i) provide the entrainment pa- 
rameters at finite temperatures ( as in Gusakov & Haensel 
(2005)) (essentially doing equation of state calculations 
"out of equilibrium" ) , and ii) develop fluid models that 
allow for the presence of both superfluid condensates and 
the quasiparticle excitations that will be present when the 
star is no longer at T = 0. 



5 Mutual friction 

A superfluid mimics bulk rotation by forming an array 
of vortices. The dynamics of these vortices, and their 
interaction with the different particle species may signif- 
icantly affect the evolution of the system. In a neutron 
star the presence of vortices leads to "mutual friction" 
between interpenetrating superfluids (neutrons and pro- 
tons). This couples the two fluids on a relatively short 
timescale, which has implications for both glitch recovery 
estimates and the damping of pulsation modes. In fact, 
an estimate by Lindblom and Mendell indicates that the 
gravitational- wave driven instability of the fundamental 
f-modes will be prohibited by the mutual friction (Lind- 
blom & Mendell 1995). The same effect is also relevant, 
albeit not quite in such a dramatic fashion, for the r- 
modes (Lindblom & MendeU 2000). 

In a recent study, we have revisited the problem of 
the mutual friction force for neutron stars (Andersson, 
Sidery & Comer 2006). It is known that the entrainment 
plays a key role in determining the strength of the mu- 
tual friction, and to make progress and study various as- 
trophysical scenarios we needed to develop a description 
of this mechanism within our multifluid framework. Our 
results provide useful checks on the classic work by Al- 
par, Langer & Sauls (1984) (who were the first to discuss 
the problem) and Mendell (1991) (whose mutual friction 
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Fig. 3 A schematic illustration of the flow around a super- 
fluid neutron vortex. Because of the entrainment effect, the 
vortex induces flow also in the protons. This leads to a mag- 
netic fleld forming on the vortex, and the scattering of elec- 
trons off of this magnetic fleld leads to a dissipative mecha- 
nism known as mutual friction. 



coefficients were used in the studies of mode-instabilities 
in superfluid neutron stars). 

Briefly, the mutual friction that is expected to be the 
most important in the outer core of a mature neutron 
star originates as follows. The superfluid neutrons form 
vortices, which represent the quantisation of the momen- 
tum circulation. That is, we have 



h 



2mn 



-J _ ijk 



(9) 



where is a normalised tangent vector to the vortex. 
Here it is important to understand that it is the circu- 
lation of momentum that is quantised, not that of the 
velocity. At first sight this distinction may seem to sug- 
gest that our model dilfers from the standard picture. 
However, this would be wrong. One must remember that 
in the orthodox description for superfluid Helium, the 
so-called superfluid "velocity" is in fact a rescaled mo- 
mentum (Prix 2004). Hence, the two pictures are con- 
sistent, the only dilference being that we make the true 
dynamical role of the variables clearer. Now, the flow of 
neutrons associated with the vortex induces flow in part 
of the protons because of the entrainment. This leads 
to the generation of a magnetic fleld, see Figure 3, the 
strength of which is estimated as 
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where pi4 = p/lO^^'g/cm'^. Finally, electrons scatter dis- 
sipatively off of this magnetic fleld, leading to a coupling 
of the two fluids in the system (the superfluid neutrons 
and a conglomerate of all the charged particles). 

To arrive at the standard expression for the mutual 
friction force, we balance the Magnus "force" , which acts 
on a vortex as the superfluid flows past it, and a resis- 
tive force due to electrons scattering off of the vortex. 
It should be noted that the Magnus effect is already 
present in the smooth averaged equations of motion (8), 
which means that in practice we only add resistivity to 
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the right-hand side of the equations. Anyway, the force- 
balance leads to the mutual friction force acting on the 

neutrons being 



/^f = BpnUv e'^''kj€kimK'^w'^p +B'pnnv e 



ijk 



(11) 



projccts±/^^ 



This force has been used to model a number of astro- 
physical scenarios. In particular, the coupling of the two 
fluids following a glitch event, see the following section, 
and the damping of oscillations in a superfluid neutron 
star. It is, however, important to understand that this is 
only a first approximation to the real mutual friction 
interaction. There are presently two alternatives that 
should be taken seriously. The first, which has been ad- 
vocated by Sedrakian and collaborators for over 20 years 
[see Sedrakian & Sedrakian (1997) and references given 
therein], is based on the notion of vortex clusters. The 
idea is that each neutron vortex is associated with a 
large collection of proton fluxtubes. This increases the 
scattering cross section for the electrons, and as a result 
the mutual friction force is many orders of magnitude 
stronger than given above. The upshot of this is (some- 
what counter-intuitively) that the coupling between the 
fluids take place on a much longer timescale (months 
rather than seconds). The other alternative accounts for 
the curvature of the vortices, and the possibility that 
they will get tangled up in a state that in many ways 
is reminiscent of turbulence. The lack of a prefered di- 
rection of the vortex array obviously affects the form of 
the mutual friction force. This is well-known from work 
on superfluid Helium (Donnelly 1991). Nevertheless, this 
alternative form for the mutual friction has only recently 
been discussed for neutron stars. Peralta et al (2005) are 
leading the way by carrying out numerical simulations 
of a turbulent two-fluid system. This work has a number 
of interesting implications, and it will be exciting to see 
where this will take us. 

Future work needs to investigate the alternative mani- 
festations for the mutual friction in detail. In my opinion 
the most pressing issue concerns the turbulence. Since it 
is known to be the key effect in experiments on super- 
fluid Helium one has every reason to expect that it will 
be relevant for neutron stars as well. 



6 DynEtmical coupling 

One of the most important problems that can be ad- 
dressed once we know the form of the mutual friction 
force regards the relaxation timescale following a pulsar 
glitch. The idea is that in the initial state the vortices are 
pinned. This leads to a build-up of a rotational lag be- 
tween the two fluids as the charged component (the crust 
and the core protons/electrons) spins down. At some 
critical relative rotation rate, the pinning force will no 
longer be able to hold the vortices in place. Global unpin- 
ning then leads to evolution described by the equations 



that we have derived (provided that the vortices remain 
straight during the evolution, a debatable assumption). 
All we need in order to be able to estimate the relaxation 
timescale following the glitch are the coefficients in (11). 
As we have recently shown (Andersson, Sidery & Comer 
2006), the dimensionless coefficients are B' = B^ and 
7^ 
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where drup = rUp — m* . Working out the dynamical cou- 
pling timescale we easily find 

n^dtp't + --- = ff 
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That is, the timescale on which the two fluids are (lo- 
cally) coupled can be estimated as 
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This suggests that the two core fluids are coupled on 
a timescale of about 10 rotation periods (for typical pa- 
rameters). This estimate should be compared to the clas- 
sic result of Alpar & Sauls (1988) who suggest that the 
coupling timescale is 400 — 10^ periods. This relatively 
rapid coupling is usually taken as evidence that glitches 
must originate in the crust superfluid, not in the core. 
Our calculations, with a coupling that is about 1 order 
of magnitude faster than the old result, obviously sup- 
port this conclusion. However, we need to understand 
why our estimate differs from that of Alpar and Sauls. 
The result may seem surprising since the parameters in 
the mutual friction force are similar in the two stud- 
ies. The difference comes from the set-up of the cou- 
pling timescale problem. Our estimate followed imme- 
diately from the two-fluid equations, and represents the 
timescale on which mutual friction couples the two com- 
ponents in the system. In contrast, Alpar and Sauls con- 
sider the equation of motion for the vortices. Setting the 
system into relative rotation, they work out the timescale 
for the vortices to relax to a new equilibrium position. 
The key difference between the two estimates is that in 
the Alpar-Sauls scenario the rotational lag between the 
fluids is maintained while the vortices relax. This means 
that their model does not conserve angular momentum 
during the relaxation, which explains why the two results 
are different. 



7 Assignments 

At the end of my presentation, I gave the audience a 
number of "homework assignments" . These represent prob- 
lems that I feel need more attention, and which I hope 



Modelling the dynamics of superfluid neutron stars 



to consider in the near future. It seems appropriate to 
conclude this brief write-up by outlining these problems. 

As we are beginning to understand the basic dynam- 
ics of the two-fluid model for neutron stars, we should 
look ahead and develop a framework that would allow 
us to make our models more realistic. Most importantly, 
we need to account for finite temperature effects and vis- 
cosity. This is a difficult problem area, but I feel that we 
cannot shirk away from it. Recently, Greg Comer and I 
have developed a flux-conservative formalism for multi- 
fluid systems, including both dissipation and cntrain- 
ment (Andersson & Comer 2006a). Our formulation has 
much in common with work on the mechanics of mixtures 
and multiphase flows. Our first study of the problem is 
quite formal, and should be relevant for any multifluid 
situation in Newtonian physics. Of course, we want to 
apply the framework to neutron stars. To do this, we 
have considered two model systems. The results make it 
painfully clear that more realistic models will be signif- 
icantly more complex than the systems we have consid- 
ered so far. Our analysis shows that the simplest "rea- 
sonable" model for a neutron star core, starting with four 
fluids (neutrons, protons, electrons and entropy) and re- 
ducing to three degrees of freedom (superfluid neutrons, 
entropy and everything else), requires no less than 19 
dissipation coefficients. 




Fig. 4 An illustration of the different kinds of motion that 
lead to dissipation via particle scattering in a two-fluid prob- 
lem (blue and red). The "red" fluid is assumed to be at rest, 
and the "blue" fluid i) flows linearly, ii) contracts and ex- 
pands, and iii) undergoes shearing motion. 



It will be a challenge to understand the role and val- 
ues of these different coefficients. Yet it may not be too 
bad a problem, because the meaning of many the of new 
coefficients is quite clear. Consider for example dissipa- 
tion due to particle scattering. In a single fluid problem 
scattering leads to viscosity due to shearing motion in the 
flow. In a situation where a relative flow is possible, there 
are more degrees of freedom that lead to particle scat- 
tering. Consider the schematic illustration in Figure 4, 
where the two fluids are represented by blue and red par- 
ticles. There are three different ways that relative flow 
can induce scattering, and result in dissipation. Taking 
the "red" fluid to be at rest, the blue fluid can i) flow 
linearly, ii) contract and expand, or iii) undergo shearing 
motion. The viscosity coefficients associated with these 
different degrees of freedom are related and should only 
differ by "geometric factors". They should certainly be 
calculable, and might actually be known from kinetic 
theory already. 
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As described in other contributions to these proceed- 
ings, it is likely that we now have the first observations 
of neutron star oscillations. This is fantastic news, but it 
provides us with a challenge that must be met in the near 
future. We need to develop accurate theoretical models 
to put contraints on eg. the crust physics via these obser- 
vations. However, if we want to do this we must consider 
fully relativistic models. This is the only way to obtain 
results that have the precision required to make a com- 
parison with observations meaningful. In fact, we might 
even have to account for the presence of superfluid neu- 
trons in the crust, a very difficult problem. Moreover, 
the observed systems are all magnetars and one might 
expect the magnetic field to be important. It seems likely 
(at least to me) that the magnetic field will couple any 
motion in the crust to the core, thus altering the nature 
of the pulsation modes and making any analysis of "pure 
crust oscillations" less relevant (Glampedakis, Samuels- 
son & Andersson 2006). This would force us to study 
global oscillations of magnetised fully relativistic stars, 
certainly not a simple problem. 

As I have already indicated, we need to worry about 
superfluid "turbulence". I would not be at all surprised 
if studies of this problem lead to results that change 
our understanding of neutron star dynamics significantly. 
One advantage is that there has been a lot of work on 
the analogous problem for superfluid Helium (Donnelly 
1991), and we can hope to benefit from these results. 
Most importantly we need to understand whether a tur- 
bulent description is relevant for neutron stars. If it is, 
how does it manifest itself? What is the effect on, for 
example, glitch relaxation and mutual friction damping 
of neutron star oscillations? 

In addition to these problems, I can think of a number 
of issues concerning multifluid aspects of exotic phases 
like hyperons and deconfined quarks. Certainly the neu- 
tron star community is bustling with exciting ideas, but 
we have a lot to do before we truly understand the dy- 
namics of superfluid neutron stars. 

Acknowledgements Much of the work described in this ar- 
ticle was done in collaboration with Greg Comer. Without his 
contribution the results discussed here would not have been 
the same. 
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